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The particle number and energy fluctuations in the system of charged particles are studied in the 
canonical ensemble for non-zero net values of the conserved charge. In the thermodynamic limit the 
fluctuations in the canonical ensemble are different from the fluctuations in the grand canonical one. 
The system with several species of particles is considered. We calculate the quantum statistics effects 
which can be taken into account for the canonical ensemble fluctuations in the infinite volume limit. 
The fluctuations of the particle numbers in the pion-nucleon gas are considered in the canonical 
ensemble as an example of the system with two conserved charges - baryonic number and electric 
charge. 



INTRODUCTION 



O . The statistical model approach turns out to be rather successful in describing the data on the particle production 
' in relativistic nuclear collisions (see e.g. Ref. |IJ and recent review @). This stimulates further investigation of the 
I . properties of the statistical models. In particular, an applicability of various statistical ensembles is an interesting 
£f*^ ' issue. The canonical ensemble (CE) Q or even the microcanonical ensemble (MCE) J] have been used in order to 
describe the pp, pp and e + e~ collisions when a small number of secondary particles are produced. At these conditions 
the statistical systems are far away from the thermodynamic limit, so that the statistical ensembles are not equivalent, 
k* - | and the exact charge or both energy and charge conservation laws have to be taken into account. The grand canonical 
CTj ■ ensemble (GCE) formulation is valid when the system volume V tends to infinity. All statistical ensemble become 
2^ \ then thermodynamically equivalent. 

The analysis of the fluctuations is a useful tool to study the properties of the system created during high energy 
' particle and nuclear collisions (see e.g. Refs. 0, 0, 0|)- An essential part of the total fluctuations measured on the 
, event-by-event basis is expected to be the thermal ones. The particle number fluctuations have been recently studied 
CD ' in the CE and MCE and compared with those in the GCE. It has been shown that these fluctuations are different 
in various statistical ensembles in the particular case of the relativistic ideal gas with a total net charge equal zero in 

■ the Boltzmann statistics approximation. The fluctuations of negatively and positively charged particles are suppressed 
in the CE 8] in comparison to the fluctuations in the GCE. This suppression remains valid in the thermodynamic 
limit too, so that the well-known equivalence of all statistical ensembles refers to the average quantities, but does not 
apply to the fluctuations. 

In Ref. @ we have studied the CE for one particle specie and zero net value of the conserved charge. In the present 
paper we extend our consideration. In the high energy proton-proton and nucleus-nucleus collisions the created system 
has some positive values of the baryonic number and electric charge. Besides, a lot of different species of hadrons 

■ are created. We study the CE particle number fluctuations (Sees. II and III) and energy fluctuations (Sec. IV) in the 
Ci ' systems with non-zero net charge and several species of particles. As the electric charge of hadrons can be both ±1 

and ±2, we consider the CE system of single and double charged particles in SecV. The effects of Bose and Fermi 
statistics are studied in the thermodynamic limit in Sec. VI. We also calculate in Sec. VII the CE particle number 
fluctuations for the ideal pion-nucleon gas which is an example of the system with two conserved charges - baryonic 
number and electric charge. We summarize our consideration and formulate the conclusions in Sec. VIII. 



II. THE GCE AND CE PARTITION FUNCTIONS AND MEAN PARTICLE MULTIPLICITIES 

Let us start with the multi-species system of +1 and —1 charged particles. In applications of the statistical approach 
to hadron production in high energy collisions the conserved charge under consideration can be the electric charge 
and baryonic number, or strangeness and charm, which are also conserved in the strong interactions. In the case of 
the Boltzmann ideal gas (i.e. the interactions and quantum statistics effects are neglected) the partition function in 
the GCE reads: 



2 



Z g . c .eXV,T, ll) 



E ••• E • 

Ni + , ATi_=0 Afy + , JVj_=0 



N 1+ \ 



N j+ \ 



n e 



(Aj+Zj)^ 4 " (Xj-Zj) 1 * 3 ' 



N- , l 



NjJ. 



II ex P( A i+^ + A J- z i) = ex P 



2z cosh 



(1) 



where j numerates the spesies, Xj± — exp(±^/T), Zj is a single particle partition function 
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and z = Zj. The V, T and /i are respectively the system volume, temperature and chemical potential connected 
with the conserved charge Q. The gj and nij are respectively the degeneracy factors and masses for the j-th particle 
species, and K2 is the modified Hankcl function. The CE partition function is obtained by an explicit introduction 
of the charge conservation constrain, Y^j — Nj-) = Q , for each microscopic state of the system: 
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Parameters Aj+ and Aj_ in the CE © are only auxiliary parameters introduced in order to calculate the mean 
number and the fluctuations of positively and negatively charged particles. They are set to one in the final formulas. 
In Eq. the integral representations of the 5-Kronecker symbol and the modified Bessel function were used [To|: 
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The averages of Nj + and Nj- in both the GCE and CE can be presented as (in the final expressions one should 
put Xj± = exp(±/i/T) and Xj± = 1 for the GCE and CE, respectively): 
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for the GCE and CE, respectively. The average number of -/V + and iV_ are equal to: 
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The mean net charge in the GCE is equal to: 

Q = (N+) g . c . e . - (JV_) 9 .c. e . = 2sinh(^) 
which leads to a simple relation which connects the values of Q and fi 
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so that 



(N±) g . c . e . = z (y + y/l + y 2 ) 
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where y = Q/2z = sinh(yLt/T). 

In the CE an exact charge conservation is imposed on each microscopic state, so that it is evidently fulfilled also 
for the average values: 
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as indeed can be easily seen from the identity I n -i{x) — I n +i(x) 
The ratios of (N±) calculated in the CE and in the GCE, 
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are shown in Fig. 1 for Q = and Q = 2. There is the strong canonical suppression effect, (N±) c . e . <C (N±) g . c . e ., for 
small systems (z -C 1), and the canonical and grand canonical ensembles become equivalent, (N±) c . e , = (N±) g _ c _ e _, in 
the thermodynamic limit z — > oo. One can see that the CE suppression effect is reduced for a non-zero net charge 
of the system as compared to a system with zero net charge. In Fig. 2 the ratios l|12|) as functions of Q — 1, 2, ... are 
shown at fixed positive values of y = Q/2z which correspond to the fixed positive net charge number densities (Q = 
corresponds to y = and is presented in Fig. 1). Small values of y mean large z, e.g. for y — 0.1 shown in Fig. 2 one 
finds 'large' z — 5 at Q = 1, so that the system is already close to the thermodynamic limit. Due to this the canonical 
suppression is small and it is the same for positive and negative particles. The case of large y differs, e.g. for y = 2 
shown in Fig. 2 the values of z are 'small' at small Q: z = 0.25 at Q = 1. The canonical suppression effect becomes 
strong for negative particles at small Q. However, the canonical suppression at large y is negligible for the average 
value of positive particle number as it should be approximately equal to Q. 
For small systems (z<l) using the series expansion 10] 
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In the large volume limit (V — > oo corresponds also to z — > oo) the mean quantities in the CE and GCE are equal. 
This result is referred to as an equivalence of the canonical and grand canonical ensembles. Using the uniform limit 
of the modified Bessel function |ltj| 
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one can easily find (note that fixed Q at z — > oo means a zero value of the net charge density and y 

(N ± ) c . e . ~ z(y+ ^l + y 2 )^ = (N ± ) g . c . e . . 
The total multiplicity of charged particles is defined as N c h 
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FIG. 1: The ratios of the mean particle numbers in the CE to those in the GCE as functions of z for Q = and Q = 2. 
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FIG. 2: The ratios of the mean particle numbers in the CE to those in the GCE as functions of Q = 1, 2, 3, ... for fixed values 
of y = Q/2z. 



III. THE SCALED VARIANCE 



A useful measure of the fluctuations of any variable X is the ratio of its variance V(X) = (X 2 ) — (X) 2 to its 



mean value (X), referred here as the scaled variance: 
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Note, that uj x = 1 for the Poisson distribution. In order to study the fluctuations of charged particle numbers the 
second moments of the multiplicity distribution have to be calculated. One finds: 
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in the GCE and CE, respectively. The scaled variances ufl ± and ch are equal to: 
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The Eqs. H26l27fl describe the particle number fluctuations of a given specie j . One can establish the general rule 
how to calculate the fluctuations of N± = J2j Nj± and N c h = N + + N-. To do this one should set Ai± = A2± = • • • = 
A± in Eqs. ill I- ill and differentiate with respect to A± in Eqs. I|5I23|I in order to get (N±) (n = 1,2). This eventually 
results in a substitution of z instead of Zj in all final formulas for the averages and fluctuations. One obtains: 
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The scaled variances uuf e and wf e calculated with Eqs. i(2T)|) and (|3"U|l are shown in Fig.0 for Q = 0, Q = 2 and in 
Fig.0]for fixed positive values of y. Using the asymptotic behavior of the modified Bessel function for z — > 0, Eq. I|13fl . 
and z, Q — * oo with y — Q/2z = const, Eqs. (|17I18|I . the limits of the scaled variances can be easily found, both for a 
given particle specie j and for the sum over all particle species: 
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FIG. 3: The scaled variances ujt.e. and u^. 130t as functions of z for fixed values of the conserved charge Q. 
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FIG. 4: The scaled variances u^e. 12!» and u£^. I30H as functions of Q = 1, 2, 3, ... for fixed values of y = Q/2z. 
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2). A large system limit z — > oo gives for fixed Q (note again that fixed Q in the thermodynamic limit z 
means a zero value of the net charge density and leads, therefore, to y = 0) 
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As one sees from Eqs. 3-4 the scaled variances reach very fast their asymptotic values. In Fig.[3]the scaled variances 
for Q = and Q — 2 can be compared (for Q = see details in Ref. 8]). One notices that their values at z — ► oo are 
the same, but the behavior at small values of z is different. Namely, if Q > 1 the fluctuations of positively charged 
particles are very small at small z, while the fluctuations of the negatively charged particles have the Poisson width. 
This can be easy understood as for small volumes the average number of positive particles is approximately equal to 
Q (see Eq. I)15|l) and the fluctuations of N+ are small. On the other hand, at small z and fixed Q the average number 
of negatively charged particles is much smaller than Q (see Eq. (IlKf) ) and the fluctuations of 7V_ are not affected by 
the conservation law. Similar physical reasons explain the behavior of the fluctuations at non-zero charge density in 
the thermodynamic limit. The Fig. 4 demonstrates the following features of the asymptotic values of uj^ e and uj~ e 
at Q ^> 1. When the charge density becomes larger (y increases) the w+ e decreases and tends to at y — > oo, while 
the uj~ e increases and tends to f at y — > oo. The physical reasons of this are seen from Eq. (|19|l which at y 3> 1 gives: 
(iV_|_) c . e . ~ z ■ 2y = Q and (A r _) c e . ~ z ■ (2y) _1 = Q ■ (4y 2 ) -1 <C Q. Therefore, at y 1 an exact charge conservation 
keeps N + close to its average value Q and makes the fluctuations of N + in the CE to be small. Under the same 
conditions, (N-) c . e . is much smaller than Q, so that the fluctuations of N- are not affected by the CE suppression 
effects and has the Poisson form. These features of the CE are in a striking difference with those in the GCE. The 
GCE scaled variances (|28ll are equal to f for 7V_, N + and N c h, and this remains valid for all values of the system net 
charge or net charge density. 



IV. ENERGY FLUCTUATIONS 



The partition function in the GCE and CE is equal to Z = J] exp(—E/T), where the sum over microstates includes 
the summation (integration) over particle momenta and summation over number of particles and over different particle 
species. Each microscopic state has the weight factor ■ exp[(fiNj+ — /iA^_)/T] in the GCE Qand S[Q — ~ 
Nj-)] in the CE J3J|. In order to calculate the average energy and its fluctuations it is convenient to rewrite the 
partition function as Z = X) ex PEj (Af+-^j+ + @j-Ej-)/T], where (3j + and /3j- are the auxiliary parameters and 
(5j + = /3j- = j3 = 1/T in the final formulas. It then follows: 
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where Zj — dZj/d/3, z^ = d 2 Zj/d/3 2 , and Zj,a±,b± are given by Eqs. (|2I6I25J) . respectively. In Eq. (f^T)f> we have intro- 
duced the average value of one-particle energy (sj) = —Zj/zj. Introducing also (e 2 ) = Zj / Zj the energy fluctuations 
can be then presented as: 
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where w J is given by Eq. (|26|l . Introducing the total energies E± = Ej± and E c h = ^2j(Ej + + EjJ), one finds: 
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fluctuations consist of two terms 



Y^{s?)zj/z. The Eas. (|44l45l) are valid in both the GCE and CE. The energy 
The first term takes into account the fluctuations of one-particle energies, the 
second one - the fluctuations of the number of particles. Most often the fluctuations of the number of particles are 
relatively more important than the fluctuations of one-particle energies. Indeed, the maximal value of the first term 
in the right hand side of Eqs. I|44I45|I is equal to (e)/3 for the particles with ra/T — > 0, and it decreases and goes to 
zero at m/T — ► oo. On the other hand, the second term in the right hand side of Eqs. (|44l45ll is equal to (e) for any 
system in the GCE. The value of ((e 2 ) — (e) 2 )/(e) in Eqs. I|44I45|I is the same for "+"and "-" particles, and in both 
the GCE and CE. The values of uj's are however different in the GCE and CE. Besides, the w+ e , u)~ e and u>° e are 
different from each other for the non-zero net charge Q. Therefore, the scaled variances of the energy fluctuations are 
different in the GCE and CE, and in the CE the values of W + , W~ and W ch differ from each other and depend on 
the net charge of the system. An example of the energy fluctuations W + , W~ and W ch for the ideal pion gas with 
Q = and Q = 2 is presented in Fig.0 One sees that the dependences of the energy fluctuations W + , W~ and W ch 
on z in the CE resemble those for ui 
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FIG. 5: The CE energy fluctuation W + , W~ and W ch in the ideal pion gas at temperature T = 120 MeV. 



The upper horizontal dotted line in Fig. 5 shows the value of (e 2 ) /(e) which corresponds to the W + = W~ = W ch 
in the GCE. The lower horizontal dotted line in Fig. 5 shows the value of ((e 2 ) — (e) 2 )/(e) which is the fluctuations 
of one-pion energy. 



V. SINGLE AND DOUBLE CHARGED PARTICLES 



In the following sections we consider the extension of the CE formalism. First, let us study the system of particles 
and antiparticles with charges ±1 and ±2. The CE partition function reads: 
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where we have used the relation exp \x (t + |)] = J2T=-oo t Ik(^ x ) . The z and z in Eq. (|46() are the one-particle 
partition functions for charges ±1 and ±2, respectively. In terms of variables c m ± (m = 1,2,4) 
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To illustrate the specific features of the considered system we present the CE results in the case Q = . As always, 
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Eqs. (|48H52|) for finite values of z and z one can effectively use Eq. (|47f) . In the thermodynamic limit V — > oo both 
z — > oo and z — > oo. In this case it is convenient to return to the integral over <f> representation in Eq. (|46|l and use it 
also for the derivatives of the CE partition function with respect to A± and A±. Using the saddle point method to 
calculate the ^-integrals one finds then for z,z^S> 1: 
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From the above formulas one finds that ~ 0.5 if z/z is either much smaller or much larger than 1. The scaled 
variance has a maximum at z = 2z . At this point one finds a->^ e = 5/9, u/^ c e = 5/6 and uf^ c e = 2/3. 



VI. QUANTUM STATISTICS EFFECTS 



It is instructive to apply a different technique to calculate the fluctuations of the thermodynamical quantities 
with the exact conservation laws imposed. This method allows to find the values of the CE fluctuations in the 
thermodynamic limit and include the effects of quantum statistics. 

The ideal quantum gas of the identical Bose or Fermi particles and antiparticles can be characterized by the 
ation numbers of the one-particle states labeled by momenta p. The GCE average values and fluctuations 
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where 7 is equal to +1 and —1 for Bose and Fermi statistics, respectively (7 = corresponds to the Boltzmann 
approximation). These expressions are microscopic in a sense that they describe the average values and fluctuations 
of a single mode with momentum p. However, the average values of all macroscopic quantities of the system can 
be determined through the average occupation numbers of these single modes. The fluctuations of the macroscopic 



observables can be written in terms of the microscopic correlator (AnJ 
correlator can be presented as: 



Afijlj.c.j., where a, (3 are + and(or) — . This 



(An«Anf) 9 . c . e . 



,,a2 



S P k S af 3 , (57) 



10 



due to the statistical independence of different quantum levels and different charge states in the GCE. The variances 
of the total number of (negatively) positively charged particles N a — E p rip are equal to: 



(AN ( 



a)g.c.e. 



a) g.c.e. 



p.k 



(K<) £ 



\ n p lg-c.e. \n k ) g.c.e. 



p.k 



We have assumed above that the quantum p-levels are non-degenerate. In fact each level should be further specified 
by the projection of a particle spin. Thus, each p-level splits into g = 2j + 1 sub-levels. It will be assumed that the 
p-summation includes all these sub-levels too. The degeneracy factor enters explicitly when one substitutes, in the 
thermodynamic limit, the summation over discrete levels by the integration: 



E 

p 



gv_ 

2tt 2 



p dp 



The scaled variance lo 

j2 



g.c.e. 

(N' al g.c 



in the thermodynamical limit V - 
- (N a ) 2 g .c.e. E P , fe <An£An£) 



reads: 



E p '< 



g.c.e. 



pi g.c.e. 



Jo P 2d P 



(58) 



'p / g.c.e. 



The Eq. H58|) corresponds to the particle number fluctuations in the GCE. To illustrate the role of quantum statistics 
let us consider the case of n = 0, i.e. (Q)g. c .e. — 0, where Q = E p a q a np*. In what follows we assume q + = 1 and 
q~ = —1, however, the formulas below are valid for any values of q + and q~ = —q + . From Eqs. 1|56|) and (|58|l one 
finds w±f ° ltz = 1 in the Boltzmann limit (7 = 0), w±f° se > 1 for the Bose particles (7 = 1) and w±ff mi < 1 for the 
Fermi particles (7 = —1). The strongest quantum effects correspond to m/T — > 0: 



±Boltz 
g.c.e. 



The scaled variance Wg h c . e . 
one finds: 



±Bose 
g.c.e. 



1.368 



±Fermi 
g.c.e. 



0.912 



6C(3) ' ~ 9ce ' 9C(3) 

for all charged particles can be easily obtained from (|58|l by replacing J2 P 



(59) 



ch Boltz 
' g.c.e. 



±Boltz 
°g.c.e. 



ch Bose 
- g.c.e. 



UJ 



±Bose 
g.c.e. 



ch Fermi 
' g.c.e. 



AiFermi 
°g.c.e. 



(60) 



The formula for the microscopic correlator is modified if we impose the exact charge conservation in our equilibrated 
system. For this purpose we introduce the equilibrium probability distribution Win^) of the occupation numbers. As 
a first step we assume that each rip fluctuates independently according to the Gauss distribution law with the average 



value (npg.ee. 



I55|) and the mean square deviation v^ 2 (|56|l : 

W(n%) oc J|exp 



2vf 



(61) 



To justify this assumption (see Ref. 0) one can consider the sum of rip in small momentum volume (Ap) 3 with the 
center at p. At fixed (Ap) 3 and V — > 00 the average number of particles inside (Ap) 3 becomes large. Each particle 
configuration inside (Ap) consists of (Ap) ■ V/ (27r) 3 3> 1 statistically independent terms, each with the average value 
(rip)g.c.e. I|55(l and the scaled variance vp* 2 (|56() . From the central limit theorem it then follows that the probability 

distribution for the fluctuations inside (Ap) 3 should be Gaussian. In fact, we always convolve rip with some smooth 
function of p, so instead of writing the Gaussian distribution for the sum of rip in (Ap) 3 we can use it directly for n™ . 

The average value of the conserved charge Q = E p a q an p is regulated in the GCE by the chemical potential /j,. If 
we impose an exact charge conservation, AQ = E p a q a Arip = 0, the distribution (|61|l will be modified as: 



W{n°) oc ]Jexp 



p , (X 



(Ar 



2vf 



\p,a / J -°° p,a 



exp 



2v 



n2 



+ iX g Q An; 



(62) 



It is convenient to generalize distribution l|62|) to W(rip, A) using further the integration along imaginary axis in the 
A-plane. After completing squares one gets: 



W(nJ,A) oc ]Jcxp 



At 



\vfq a ) 



2vf 



vfq a2 



(63) 
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and the average values are now calculated as: 

fZo^f-ooTl P , a dn ? W(n2,\) 
The Eq. (|64[) has the meaning of the CE averaging in the thermodynamic limit V — > oo. One easily finds 

((An* - vf\q a )(A4 - «fV)> = S pk 6 a0 vf , (A 2 ) = - fe ' « An P ~ < A< ? Q )A) = 

\p,a / 

so that it follows: 



(64) 



(An^Anf) = S pk 5 af} vf - vfq a wfV (A 2 ) + (An™ A) vftf + (Anf A) vfq a (65) 
= <5 P fc S aP vf + vfq a vfqP (A 2 ) = 5 pk S aP vf 



v a2 n a v^nf 3 



Y v a2 q a2 



By means of Eq. (|B3|l we obtain: 



E p vf (E P < 



(TV 2 ) - (a q ) 2 i. ,^„<V 



2 



(66) 



The Eq. (|55|l leads to t>™ 2 = (n™) 9 . c . e . in the Boltzmann approximation, so that from Eq. I|66(l one finds (y = Q/2 
sinhOu/T)): 



a = 1 exp(a^/r) 1 y 

exp(/x/T) + exp(-/VT) 2 2^1^ ' 1 j 

which coincides with Eq. Formula can be obtained from after replacing ^ p — > X) P a > ana ^ ^ ^ s ^ ne 

same as Eq. Ij39(l . At // = from Eq. (|66|) we find the CE scaled variances: 

ch Boltz r> ±Boltz ch Bose r> ±Bose ch Fermi r> ±Fermi fac\\ 

U c.e. — ZuJ g.c.e. > w c.e. — ^ w c.e. ) W c.e. — ZuJ c.e. • (.^J 

As seen from Eqs. I|59I68I1 the scaled variance of (negative) positive particles with Bose or Fermi statistics in the 
CE is as half as large as the corresponding scaled variances in the GCE. Therefore, the CE suppression of the particle 
number fluctuations in the thermodynamic limit works at ji = in the quantum systems similar to that in the 
Boltzmann case. This result can be rephrased in another way: the Bose enhancement and Fermi suppression of 
the GCE fluctuations remain the same in the CE for the ui^ e at p = in the thermodynamic limit. The Eq. (|69l) 
demonstrates that the scaled variances of all charged particles in the CE for any statistics are by a factor of 2 larger 
than the corresponding scaled variances for (negative) positive particles, whereas in the GCE these scaled variances 
presented by Eq. (|6Tj|) are equal to each other. 

Comparing Eq. (|65|l and Eq. (|57|l one notices the changes of the microscopic correlator due to an exact charge 
conservation. Namely, in the CE the fluctuations of each mode is reduced, i.e. the ((An™) ) calculated from Eq. (|65l) 
is smaller than that in Eq. Ij56(l . and the anticorrelations between different modes p ^ k and the same charge states 
a — (3 appear. These two changes of the microscopic correlator result in a suppression of the CE scaled variances u)" e 
in a comparison with the GCE ones uj^ c e (compare Eq. (|66l) and Eq. (|58|) ~). Therefore, the fluctuations of both A_ 
and N + are always suppressed in the CE. As we have seen in the previous sections the behavior of N c h fluctuations in 
the CE can be more complicated. This is because of the correlations of different modes p k for the different charge 
states a = —(3 (i.e. the second term in the right hand side of Eq. is positive for a = —(3). 

The exact charge conservation should also lead to the canonical suppression of (n™), and this should result in the 
canonical suppression effects for (N a ). They are, however, absent in the present formulation, so that formula (|65|) 
for the microscopic correlator is not enough to calculate (A" 2 ) and (N a ) 2 separately with an accuracy corresponding 
to the effects of the canonical suppression. Nevertheless, it does allow us to calculate their difference ((AJV a ) ) 
with the effects of the CE correctly included. This means that the canonical suppression effects in the occupation 
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numbers (n") lead to the changes of the order of (N a ) in both (N%) and (N a ) 2 , but these changes are the same 

and the correction terms are cancelled in the calculation of ((AA Q ) 2 ). Therefore, the macroscopic fluctuations of 
multiplicities are not affected by the CE corrections to the average particle numbers. The scaled variances of the CE 
in the thermodynamic limit V — > oo feel the consequences of an exact charge conservation due to the suppression of 
the single mode fluctuations ((Afi™) ) and due to the (anticorrelations ) correlations between different modes p ^ k 
with the (same) different charge states a,/3. All these effects are absent in the GCE. 



VII. A SYSTEM WITH TWO CONSERVED CHARGES 



In the previous sections we have considered the system with one conserved charge. In high energy collisions the 
measurements of fluctuations for the particle numbers and (transverse) energies are mainly done for electrically charged 
hadrons. Therefore, in applications of the CE results to an analysis of the data on fluctuations it would be reasonable 
to start with the case when the charge Q is assumed to be an electric charge. On the other hand, other conserved 
charges are also present in the system created in high energy collisions. In this section we consider the system with two 
exactly conserved charges - electric charge Q and baryonic number B. As an example we study the ideal pion-nucleon 
gas and neglect the quantum statistics effects. This is the simplest realistic case where we can study the influence of 
an exact B conservation to the CE fluctuations of the electrically charged particles. The partition function of this 
system in the CE is: 



Z c . e XV,T,Q,B) = J2 E 



E 

JV„,JV B =0 N n ,N n =0 N^.Nfr-- 



(XpZp) P (XpZp) P (X n Z n ) " (XfiZn) 



NJ. 



N l 

ly n- 



(X^+Z n ) n+ (A^-Zt,-) 



x S [ (N p - N p + N„+ - JV) - Q ] 6[{N P - N p + N n 



2tt 



2tt 



<kp . .. 

7TI / 7TL exp( 



iB, 



x exp 



-N n )-B] 

A<P+4>) _j_ ^_ e -i{<P+<t>) 



x exp [z n (A„ e* + A fi e"**)] X exp [z^ (A T+ e iv + X^- e"^)] = ^ h- Q {2z p ) I k+B - Q (2z n ) I k (2z*) , (70) 



k=—oo 



where we have used that exp[x(£ -I- 1ft)] — Y^kLo^ k ^(2x) - From Eq. JTHJ) it follows: 



(A? ) 



where j numerates pion, neutron and proton, a = 1 corresponds to particles 7r + , n, p and a 
ir~ , n,p, and (m = 1,2) 



(71) 

T to antiparticles 



771, a 



J2 Ik+ a - m - Q {2z p ) h+B-o&Zn) 4(2z w ) x [Z c , e .(V,T,Q,B)] 1 , 



k— — c 



C, Q = E h- Q (2z p ) I k+a . m+B ^ Q (2z n ) I k {2z,) x [Z c . e .(V,T,Q,B)Y 

k=—oo 



E 

fc— — oc 



I k + a - m -Q(2zp) I k+a . m+B -Q(2z n ) I k {2z v ) x [Z c . e .(y,T,Q,B)]~ 



(72) 
(73) 
(74) 



Formulas for the cross-averages (NiNj) can be obtained in the similar manner. The calculations with Eqs. (|72I74|) 
are effective for small systems. In this case the fc-sums in the above equations converge rapidly and small number 
of terms lead to the accurate results. In the limit of large system volume we can use another technique, similar to 
that developed in the previous section. This leads to simple analytical results. Using this method one can obtain, 
for example, the scaled variances for (negatively) positively charged particles in the thermodynamic limit. The same 
pictures can be obtained directly from Eqs. I|72l74f> by numerical calculations at z pi z n , z^ 1 . 

First, we consider the case when the electric charge Q is exactly conserved and the baryonic number B conservation 
is treated within the GCE. This results in: 



w± = 1 



4 + z t 



(75) 
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where 



zf = Zj exp (±^) , x 3 = exp (^) + exp (-^) , (76) 



and the chemical potentials ^ are equal to fi p = fiQ + fis for protons and [i^ — [iq for 7r + — mesons. When both 
Q and B are exactly conserved, the CE scaled variances of (negatively) positively charged particles are equal to: 



%<n (Zji^n ~t~ ^7t-^7t) ~T" (ZpXp -\- Z n X n ) -\- 2z Z n X n 

^Q.b = 1 - f ± - ± w 7 T ~ v » ( 77 ) 

\Zp ~r 2^7r J \ZpXpZ n X n -r Zp^p-2'7r*£'7r ~r Z n X n Z^X^ J 



where x n = exp(fig /T) + exp(~/ig /T). Let us repeat that both luq l|75|) and luq g (|77|l are obtained in the thermo- 
dynamic limit V — > oo. The (|75|l corresponds to the CE for electric charge and the GCE for baryonic number. 
The luq b <|77|) corresponds the CE for both conserved charges. It is easy to prove that LUq B < loq , i.e. an additional 
exact conservation law reduces the fluctuations. However, the additional CE suppression of the scaled variances of 
(negatively) positively charged particles due to the exact baryonic number conservation is rather small. We have 
plotted lf?5|l and (fTTjl in Fig.Elfor /iQ = to study the influence of baryon charge conservation on the fluctuations of 
electrically charged particles. As one can see from Fig. 6 the exact CE baryonic charge conservation leads to a little 




FIG. 6: The scaled variances u}q b (dashed line) and lOq b (dashed-dotted line) given by Eq. 1771 . The dotted lines show the 
scaled variances ujq and ujq given by Eq. 1751 . The solid line presents the scaled variance for all charged particles Uq . The 
results correspond to T = 120 MeV. 



additional suppression and does not change the result lo + = uj~ = 0.5 for zero values of the baryonic and electric 
net charges. Moreover, one can prove that at zero net charges any ideal Boltzmann gas with two exactly conserved 
charges (i.e. for any combination of particle charges and their masses) leads to the scaled variances equal to ui^ = 0.5 
in the thermodynamic limit. 

On the other hand, Fig. 6 demonstrates a strong dependence of the ujq and ujq values on the net baryonic density, 
it is not important whether the baryonic number treated within the CE or the GCE. The matter is that in the pion- 
nucleon gas the electric charge equals to Q = N p — N p + N n + — N„-. At [1b — the electric charge of the system 
is close to zero. Then one finds lOq ~ luq ~ 0.5 (compare to Fig. 4 at y = 0.1). The hb > leads to (N p ) > (Np), 
and this means a non-zero electric charge of the system. In this case an exact electric charge conservation leads to 
ujq > luq (see Fig. 4). At [Ib ^ T the electric charge density becomes large due to (N p ) / (Np) 3> 1, so that luq — > 
and luq — > 1 (compare to Fig. 4 at y = 2). 
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VIII. SUMMARY AND CONCLUSIONS 

We have considered the particle number and energy fluctuations for different systems within the canonical ensemble 
formulation. The results are compared to those in the grand canonical ensemble. We have studied the system with 
arbitrary number of different particle species and non-zero net charge in Sees. II and III. An exact charge conservation 
reduces the values of N + and N- fluctuations in the thermodynamic limit. At the non-zero net charge Q the canonical 
ensemble predicts a difference for the fluctuations of N + and N_ , they also different from the fluctuations of all charged 
particles N c h = N + + N- . All these features of the canonical ensemble are in a striking difference with those in the 
grand canonical ensemble. We have demonstrated in Sec. IV that the energy fluctuations of the system are mainly 
determined by the fluctuations of the number of particles and have the same volume dependence. Therefore, the 
energy fluctuations are rather different in the canonical and grand canonical ensembles. We extend our canonical 
ensemble results and calculate the particle number fluctuations in the system of single and double charged particles 
in Sec. V, include the quantum statistics effects in Sec. VI, study the systems with two conserved charges in Sec. VII. 

The canonical ensemble suppression effects for the charged particle multiplicities are well known, and they are 
successfully applied to the statistical description of hadron production in high energy collisions The canonical 
ensemble formulation explains, for example, the suppression in a production of strange hadrons and antibaryons 
in small systems, i.e. when the total numbers of strange particles or antibaryons are small. This consideration 
demonstrates a difference of the canonical and grand canonical ensembles - the statistical ensembles are not equivalent 
for small systems. When the size of the system increases all average quantities in both ensembles become equal. It 
means that in the thermodynamic limit V — > oo the canonical ensemble and grand canonical ensemble are equivalent. 
Results of Ref. || and the present study demonstrate that there are also the canonical ensemble effects for the 
fluctuations. In contrast to the canonical suppression of average multiplicities, the canonical effects for the multiplicity 
fluctuations do survive at V — > oo and they are even most clearly seen in the thermodynamic limit. The changes of 
the scaled variances due to an exact charge conservation of the canonical ensemble are not small (about 50 percent 
effects) and they are in general different for positively, negatively and all charged particles. To observe these new 
canonical ensemble effects in an analysis of the data on multiparticle production, several points should be clarified. 
To use the condition of an exact charge conservation one has to apply it to the system of all secondary hadrons 
formed in high energy collisions, and this should be done on the event-by-event basis as we are interested in the 
system fluctuations. In the experimental study only a fraction of produced particles with the conserved charges is 
detected. Introducing a probability q that a single particle is accepted in the detector a simple relation between the 
scaled variance of the accepted particles, uj acc , and the scaled variance of all particles in the statistical ensemble, u, 
was obtained @: uj acc = q • w + (1 — q). To observe the real event-byevent fluctuations u one needs q ~ 1, otherwise 
if q -C 1 one always obtains u acc ~ 1 and makes a wrong conclusion that the fluctuations correspond to the Poisson 
distribution. This fact is of a very general origin, and because of relatively small experimental acceptance a large 
part of the event-by-event fluctuations in high energy multiparticle production is lost. To observe many interesting 
event-by-event fluctuations, for example, due to the QCD critical point (see, e.g. 0] an d references therein), one 
should accept an essential part of all secondary particles. In this case the role of an exact charge conservation would 
increase. It would also have a strong influence on an extraction of the so called dynamical fluctuations (see Ref. 0) 
from the event-by-event data. 
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